Localization in a rough billiard: A a model formulation 
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We consider the quantum dynamics of a particle in a weakly rough bilhard. The Floquet operator 
for reflection at the boundary is obtained as a unitary band matrix. The resulting dynamics in 
angular momentum space can be treated in the framework of the one- dimensional supersymmetric 
nonlinear a model. We flnd analytically localization and the corresponding localization length 
5 = Del where Dd is the classical diffusion constant due to boundary scattering. 
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There is a close relation between random matrix 
theoryEI and Efetov's supersymmetric nonlinear a model 
for disordered metalsB. Efetov calculated from the a 
model in the zero dimensional limit (energy scales smaller 
than the Thouless energy) the random matiis level statis- 
tics. Later Bohigas, Giannoni, and SchmitEla studied the 
level spacing statistics of chaotic billiards by numerical 
and semiclassical methods also confirming the random 
matrix behavior. Only, quite recently, Muzykantskii et 
al.El and Andreev et alB rederived the a model (contain- 
ing the zero dimensional limit) for this type of systems. 

The chaotic billiards considered in Refs. pM were 
characterized by one typical time scale t^, the time be- 
tween collisions with the boundary. There are also other 
type of chaotic systems in which a second time scale 
tn ^ tc appears, to is the time to cover the acces- 
sible phase space diffusively. A famous example is the 
kicked rotatorlZIa which exhibits diffusion in the momen- 
tum space (and localization in the quantum case). Aply.- 
ing a new technic for averages over unitary—matricesa, 
Altland and Zirnbauer were able to derivelij the one- 
dimensional (T model for this example. In this way, they 
obtained an analytical proof of localization and, for short 
time scales, also of diffusion in the kicked rotator model. 

Recently, the diffusive regime was also identified in two 
types of chaoticj-billiards: a nearly circular Bunimovich 
stadium billiardO and slightly rough Jailliards with a gen- 
eral smooth boundary perturbationt3. In these exam- 
ples, the classical diffusion constant in angular momen- 
tum space leads to an estimate of the localization length 
for the quantum case. If this length is below the classical 
angular momentum boarder at given energy, the classi- 
cally chaotic billiard is no longer ergodic in the quan- 
tum case. Then the levelpSnacing statistic changes from 
Wigner-Dyson to PoissonOllS. 

In this paper, we consider as in |13] a nearly circular 
billiard with a rough surface defined by the angle depen- 
dent radius R{9) = Ro + AR{6). The deformation is cho- 
sen as a function that depends only on a finite number of 



harmonics AR{9)/Ro 
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random complex coefficients. This type of billiard may 
have important physical applications in contrast to the 
stadium billiard which is of more theoretical interest. As 



examples, we can mention surface waves in water dcoplets 
which are practically static for liglit propagationllj, non- 
ideal surfaces in microdisk lasergl^, apd capillary waves 
on a surface of small metallic clustersE^I. 

The classical diffusive behavior due to boundary scat- 
tering is essentially determined by the roughness k{9) — 
{dR/d6)/Ro and the diffusion constant in angular mo- 
mentum space is given byli3 
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with the angle average k^ ~ {K^(0))g. l,nax — rnRoVp/h 
is the classical angular momentum border at given ve- 
locity vp and li is the initial angular momentum. The 
diffusion constant M) corresponds to an integer time vari- 
able t measuring the number of collisions, a notion which 
we also adop|t below. Using the analogy to the kicked ro- 
tator modelB and numericaLxirecise calculations of eigen- 
states, it was demonstratedE£l that the localization length 
is given by ^ = Dd- 

In this paper, we concentrate on an analytical approach 
for the quantum case which relies on certain simplifica- 
tions but is otherwise not restricted in the parameter 
range as the direct numerical method. We will first de- 
rive an expression for the time evolution or Floquet op- 
erator corresponding to the refiection of the particle at 
the boundary. In the second part of the paper, we map 
the quantum dynamics of this operator onto the super- 
symmetric nonlinear a model. Here we will use a new 
technique which is quite different but presumably equiv- 
alent to Zirnbauer's methoda. 

To investigate the quantum problem, we expand the 
wave function ?/'(r, 9) with energy E — h k'^ /2m in terms 
of the Hankel functions which form a complete set: 



^{r,9)=Y,{aiHl^\kr)e 
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hiH\-\kr)e'' 
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In principle, the regularity of -0 at r = requires ai — bi 
and the coefficients ai are therefore the amplitudes in 
angular momentum space. However, for conceptual rea- 
sons we disregard this condition and concentrate on the 
boundary condition 'ip[R{9),9] = 0. This condition re- 
sults in a second equation bi = J2i' ^u' i^) '^i' where 
Sw {E) has the meaning of a scattering matrix for a wave 



packet reflected at the rough boundary. For time scales 
smaller than the collision time the particle does not see 
the origin ?■ = and the dynamics are determined by 
an effective scattering problem with the "Floquet oper- 
ator" S{E). The energy eigenvalues of the full billiard 
are determined by the condition detfl — S{E)] = 0. To 
calculate S{E), we multiply Eq. (||) at r = R{9) by 

r I 

obtain S = -(Af"')"! A(+) where A^^' (s == ±) are ma- 
trices given by 



Furthermore, we can consider the "quantum" diffusion 
constant defined by Dq — J2ii^ ~ 0^ I'^FzP- ^^ ^^'^ limit 
Unax — > oo, Z?q and Del coincide due to the completeness 
of the angular momentum states |Z >. 

The phases /i/ are rapidly varying with / and in the 
regime pi classical chaos they can be considered as quasi 
randonJlH. We are thus led to the effective random matrix 



and integrate over 0. From this we model S;~; = e'^''+'^" Uo(l-l) with independent random 
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dO €*('-')'' hI+^ [kR{0)] H^f [kR{0)] . (3) 
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In the following, we will consider small values of the 
roughness, k <C 1, and therefore small boundary per- 
turbations AR <C Rq. We also concentrate on the semi- 
classical limit 1 <C k AR ^ k Rq, corresponding to a 
wavelength smaller than the geometrical length scales of 
the billiard. In this limit, there is an effective cutoff in 
the angular momentum space |^| < Imax ~ kR^. The 
modes |/| > Imax correspond to evanescent modes which 
only give weak contributions. For the propagating modes 
I'l < Imaxi we can use the quasiclassical approximation 
of the Hankel functions: 

H\^\kr) w 2 [2TTki{r)r]-^/'^ exp[±i(/Xi(r) - 7r/4)] (4) 

where ki{r) ~ fc(l — rf/r'^Y/'^, fJ,i{r) ~ J^ dfkiir) = 
ki{r)r—\l\ arctan[fci(r) r/|Z|] and n « |Z|/fc is the classical 
turning point. This approximation is excellent for all r ^ 
ri in contrast to the standard asymptotical expression 
(which is only correct for r S> n). The integral in Eq. (g) 
can in principle be evaluated by a saddlepoint approxima- 
tion. For the case s = — , we immediately obtain the diag- 
onal approximation Aj' « 5ii 4/[fc;(i?o)^o] whereas for 

s = -I- there is for each value oil — I a set of saddlepoints 
Os determined by l~l = {ki[R{e s)]+ki[R{e s)]} AR'lOs) « 
s) = \/Dci {k{9s)/k) if both I and I are 
This expression is just the classical 
I #,ven by the ideal reflection 
For \l -l\ > \/I5d there 
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close to some value h 

map for one collision I 

at the boundary at angle 9J^. 

are no saddlepoints Og and tiie matrix A-'^' (and thus S) 

is an effective band matrixEJ of width ^^ y/Dd- For the 

subsequent calculation, the semiclassical expression of S 

in terms of the 9s is not very convenient, and we keep 

the 9 integral for A^+^ instead. Expanding the phases fii 

for small A_R and choosing values of / and I close to k, 

we finally obtain 

g^ r^ ^ifiliRo) + ifJ.i{Ra)+i-rr/2 _^ h^i2 ki^{Ro) AR{e)u ^ /-g-j 

This expression provides the first key result of this paper. 
Obviously, the matrix 5' is unitary with a form similar 
to the Floquet operator of the kicked rotator if one re- 
places Hi ^ tP and 2 ki-^Rg) AR{9) -^ /ckr cos{9) where 
T and fckr are the standard kicked rotator parameters. 
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tion m = Re j:z= 



--< l\ exp[i/(6')] |; >. Thefimc- 
gime jg ^ general kick potential 



with a finite number of harmonics [am = '^\/lmax ~ ^i Im 
for the rough billiard case (^]. For later use, we intro- 
duce the notations p — {f{9)'^)e — 1/2 Yl,m \^m\ and 
D = {f{9f)e = 1/2 Y.m ^"^ l"ml for the typical value 
of / and the "diffusion" constant D. In the following, 
we focus on the case /^ ^ 1 and consider the average 
probability W{l,lQ]t) = {\ <l\S* %> \^) for a transi- 
tion Zo ^^ I after a time t. The subscript /i denotes the 
average over the random phases. The Laplace transform 
of W is given by 

= — r dE(< ^oKl + e-*(^-"/2-*o) ^t)-i |; > 
27r J-TT \ 

X <;|(l + e'(^+^/2+'")5)-i|;o>\ . (6) 

The average over E can be absorbed in the phase av- 
erage and therefore we put E ~ 0. In order to apply 
the supersymmetric technique, we write W{l,lo;uj) = 
dzdz* F{z,z*)\^^^,^^ with the generating function: 

F(z,z*) = /sdet-i[l-fe*("/2+'°)C/ + J(z,z*)]\ , (7) 
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J{z) 



J{z) 



(8) 



J(z) = zEu^ ® Pb, U = e^^^C/o = e^'^S'e-^^. Here 
we have embedded the matrices in a super matrix space 
with two additional gradings, advanced-retarded grading 
(AR) and fermion-boson grading (BF). The blockstruc- 
ture in (||) corresponds to the AR-grading. Ei^i^^ is a ma- 
trix in ^-space with only one nonvanishing entry 1 at the 
position {l,lo). Pb is the projector on the bosonic sub- 
space and the notation for U , Uq, S, e^^ refer to operators 
in the I space. In contrast to conventional applications 
of the supersymmetry method, we are not dealing with 
Gaussian variables but withjtandom phases instead. Re- 
cently, Zirnbauer introducedl3 for this case a new type of 
Hubbard-Stratonovich transformation which relics on in- 
teresting and also quite involved mathematical concepts. 
Here, we proceed differently by replacing the random 
phases in a suitable way with gaussian variables: 



^2i{fii+iQ) 



1 - 2iAl {Hi + iO + iAi Aj)-^ Ai (9) 



where iJ; is an iV x TV gaussian random matrix (of the 
orthogonal ensemble in the limit iV — > cxd) with vari- 
ances {H'^~) = (l + (5,, ~^/N. Ai [A]) is an iV-dimensional 

column (row) vector such that A\ Ai ~ \ and Ai A\ is 
the projector on one particular state (e.g. j = 1) in 
the Hilbertspace belonging to Hi . The expression on the 
r.h.s. of ^ describes the 1 x 1-scattering matrix ot-a 
chaotic cavity whichisidcally coupled to one channcllij. 
It is well establishedoo that the distribution of it (for 
the above given values of Ai and (i?^r)) is given by the 
circurlar orthogonal ensemble which reduces to a uni- 
formly distributed random phase in the 1 x 1-case. The 
subscript I denotes the fact, that the states | j, I > belong- 
ing to each Hi are by definition completely independent 
(orthogonal) for different values of I, i.e. the Hi act on 
different l-suh blocks. In the following, we denote by H 
the operator on the full set of states for all I and j which 
contains the Hi in the diagonal l-blocks. In a similar 
way, we can arrange the operators Ai, {Aj) to operators 
A {A') which perform transitions from the J7-space to 
the iJ-space (iJ-space to C/-space). 

The generating function can be expressed as a superde- 
terminant in the iJ-spacc: 

F{z,z*) = (sdet-f^i\^[H + iOA + AB{to; z, z*)A'']\ 

xSdetli^{l + Uo + J) , (10) 

' B{uj,z) 

B{-Lu,zy 

B{uj, z) = i[l - e"^/"^ Uo + Jiz)] [1 + e"^/2 ^^ ^ J{z)]-\ 



B{uj; z, z*) = 



(11) 



Here the subscripts (Z,j) or (I) indicate the type of space 
over which the superdeterminant is taken. A = CT3 
is the third Pauli matrix in the AR-grading. The en- 
semble average is now performed with respect to the Hi . 
In order to understand the physical meaning of this re- 
sult, it is instructive to consider the matrix B at van- 
ishing source term and frequency: Bq — -6(0; 0, 0) — 
i(l — Uq)/{1 + Uq) = tan[/(6')/2] which is hermitian and 
diagonal in the ^-representation. The generating func- 
tion ( [l0|) corresponds to an iV-orbital Hamiltonian with 
blockdiagonal diagonal GOE-entries and a coupling be- 
tween the blocks given by ABqA\ For the moment, we 
consider the case of orthogonal symmetry Bj = Bq cor- 
responding to f{9) = f{—0). This case is actually not 
generic for general rough billiards. Below, we discuss the 
more generic unitary case. 

The technique to derive the supersymmetric riojilinear 
a model for this type of systems is well knownE3. Per- 
forming the standard steps, we obtain: 

F{z,z*) = Sdet^^il + Uo + J) /i^Qe-^IQl, (12) 



C[Q] = -Stm^HBico; z, z*) + iQ] 



(13) 



entries Q{1) that are 8 x 8-supermatrices with the non- 
linear constraint (5(?)^ = 1. The Q{1) are taken from 
the space that, corresponds to the a model with orthogo- 
nal symmetryo'Ej. Omitting the details of the (standard) 
derivation of (K), we mention however two points: First, 
the particular form of the operators A and At allowed us 
to perform the trace over the j-quantum number on the 
level of the a model action. In this trace, for each I 
only the particular state with ) = 1 contributes result- 
ing in the action (|3|). Second, the necessary saddlepoint 
approximation to arrive at the nonlinear constraint be- 
comes indeed exact due to the limit iV — *■ cxd. 

We are interested in the limit of long times (i ::^ 1) 
or correspondingly of small frequencies {to <^ 1). Fur- 
thermore, the coupling term Bq favors slow fluctuations 
of the Q{1) with I. We have therefore to apply a gradi- 
ent expansion on the action with the first nonvanishing 
contribution: 



1, 



1 



Str(,) ln(Bo + *Q) ^ gStr(,) (((1 + B^)-! [Bo, Q]) 



We keep the linear order in oj and use [g{6),Q{l)] ~ 
i g'{6) diQ{l) for slowly fluctuating Q{1) and some func- 
tion g{9). Going over to the continuous limit in l- 
space, we finally obtain the standard action for the one- 
dimensional a model 



C[Q] 



1 
32 
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Ai{u; + iO)StT{QA) 

(14) 



Here Q is a block diagonal supermatrix (in /-space) with 



+ source term contributions where D = {[f'{0)]'^)e{= 
Dq = Del) is the above introduced diffusion constant. 
Concerning the source term, we mention that the Sdet- 
prefactor in (n2|) requires some care. A priori, it generates 
Greens functions of the type {1 + Uq)^^ which correspond 
to the "clean" system without disorder and do not show 
any localization. However, a careful analysis shows that 
this contribution cancels with a corresponding term from 
the action. 

For the unitary case, we have to decompose Bq in 
even and odd contributions in 9 resulting in Bq = 
-Bo,ovcn + TaSo.odd with T3 = aJ in the T-grading for 
the time reversel symmetry (this grading appears in ad- 
dition to the AR- and i?i^-grading when deriving the 
CT- model) . Assuming that the odd term is comparable to 
the even term, we obtain in the action the extra term 
~ J dl Sir {[Q^Ts]"^). For slowly fluctuating Q{1), this 
term produces massive modes in the offdiagonal T-blocks 
of Q. These "cooperon" modes ajre therefore suppressed 
and Q becomes T-blpck diagonalcl. As a result, we obtain 
the unitary a modeS with the same action (|lj) . 

It is well known that the a model (|lj), gives rise to 
diffusive dynamics with the diffusion constant D for suf- 
ficiently short length scales. For long length scales it 
provides as in quasi Id disordered wiresESH localization 
with the localization length ^ = PD/2 where /? = 1 (2) 



is the symmetry index for orthogonal (unitary) symme- 
try. The locahzation length for /? = 2 agrees with the 
numerical findings of Ref. |Q . 

The above used gradient expansion is well justified for 
the limit /^ 3> 1. It is an interesting question in how far 
the action (Q) remains correct for p < 1. The answer 
is indeed not obvious, since we may still have D ^ \ for 
A/ 3> 1 . Now we can apply on ( |l3| ) an expansion for small 
Bq ~ (1/2) f{9). The "kinetic" part of the action is then 
given by a term ~ ^^ j l"^u-f| P Str[(5(Z) Q{1)]. This con- 
tribution couples the Q-fields on a strip of width M . The 
excitation "energy" due to Q{1) fluctuations with a wave- 
length M (the "smallest" nonzero mode) is estimated as 
~ M{D/M^) ~ D/M. The gradient expansion is justi- 
fied if this mode gives only a small contribution, i.e. if 
it becomes massive for D :^ M. For D <^ M, the Q{1) 
fields may fluctuate nearly independently and we enter 
a regime in which the eigenphases and eigenfunctions of 
the matrix S can be calcuated by simple perturbation 
theory due to weak coupling. 

In [12], a model with harmonics of different weights 
ctm ~ 1/m was considered such that D ~ M/^. Here, 
the two conditions /^ > 1 and D > M coincide and we 
have only two regimes of weak and strong coupling b^ 
tween different Z-states. However, we may also stud}^ 
a model with harmonics of equal weight, a^ ~ const. . 
Then we have D ^ M"^ p and there is a third nontrivial 
regime charaterized hy M < D < M^ where the a model 
action (|lj) is valid and the typical coupling is still small. 
This regime correponds for 5 to a band matrix w^ 
strong diagonal in the so-called Breit-Wigner regime 

In summary, we have derived the Floquet operator in 
angular momentum representation for a weakly rough bil- 
liard. This operator describes the reflection of a quantum 
particle at the boundary. In the chaotic regime, the dy- 
namics can be mapped on a supersymmetric nonlinear a 
model in the one dimensional angular momentum space. 
As a result, we obtain analytically diffusion with the dif- 
fusion constant Dd [Eq. (|^)] for short time scales and 
localization with the localization length ^ — Dd- A gen- 
eralization to three dimensional billiards is certainly an 
important and interesting topic of future research. 

We have also introduced an alternative way to Zirn- 
bauer to deal with the random phases. This technique 
relies on well known relations between the circular and 
gaussian random matrix ensembles and can also be ap- 
plied to random unitary matrices of arbitrary dimension. 
In analogy to Altland and Zirnbauer, it is also possi- 
ble only to average over the quasi energy as indicated 
in (|^). In this case, the matrix Q in (|2|) will be non- 
local in I and the subsequent expansion of the action 
will be more involvedt2l. However, we do not expect 
any serious modiflcations from this. Especially, we re- 
cover for f{d) ~ fckr cos{9) the asymptotic expression 



e = klJ4 fo, 



localization length of the quantum 



kicked rotator 
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